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ABSTRACT 

Context. Spectral observations of classical T Tauri stars show a wide range of line profiles, many of which reveal signs of matter inflow 

and outflow. Ha is the most commonly observed line profile owing to its intensity, and it is highly dependent on the characteristics of 

the surrounding environment of these stars. 

Aims. Our aim is to analyze how the Ho- line profile is affected by the various parameters of our model, which contains both the 

magnetospheric and disk wind contributions to the Ho- flux. 

Methods. We used a dipolar axisymmetric stellar magnetic field to model the stellar magnetosphere, and a modified Blandford & 

Payne model was used in our disk wind region. A three-level atom with continuum was used to calculate the required hydrogen 

level populations. We used the Sobolev approximation and a ray-by-ray method to calculate the integrated line profile. Through an 

extensive study of the model parameter space, we investigated the contribution of many of the model parameters on the calculated 

line profiles. 

Results. Our results show that the Ho- line is strongly dependent on the densities and temperatures inside the magnetosphere and the 

disk wind region. The bulk of the flux comes most of the time from the magnetospheric component for standard classical T Tauri star 

parameters, but the disk wind contribution becomes more important as the mass accretion rate, the temperatures, and the densities 

inside the disk wind increase. We also found that most of the disk wind contribution to the Ho- line is emitted at the innermost region 

of the disk wind. 

Conclusions. Models that take into consideration both inflow and outflow of matter are a necessity to fully understand and describe 

classical T Tauri stars. 

Key words. Accretion, accretion disks - Line: profiles - Magnetohydrodynamics (MHD) - Radiative transfer - Stars: formation - 
Stars: pre-main sequence 

1. Introduction CTTS yielded strengths on their surfaces on the order of 



10 3 G persistently over time scales of years (Johns-Krullet al. 

Classical T Tauri stars (CTTS) are young low-mass stars 199 9; SvmingtoHSlQOS; Johns-Kmll 200j). These fields are 

(< ZMe), with typical spectral types between F and I M, which h tQ di ^ circumstellar disk at some stellar 

are still accreting material from a circumstellar disk. These stars radii QUtside ±e steUar surface classical T Tauri stars also ex _ 

show emission lines ranging from the X-ray to the IR part of ^ inyerse p Q { (IpQ ffles> mosfl - m ±e Balmer ^ 

the spectrum, and the most commonly observed one is the Ha ^ ^ from ^ w her excitation leveK and some metal lines 

. £ line. Today, the most accepted scenario used to explain the ac- whkh are directl Hnked tQ the MA phenomenon (Edwards et al. 

V , cretion phenomenon in these young stellar objects is the magne- 1994) These ^ m broadened to velocities on the order of 

tosphenc accretion (MA) mechanism. In this scenario, a strong hundreds of kms -i, which is evidence of free-falling gas from a 

Ctf ; stellar magnetic field truncates the circumstellar disk near the distance of some stellar radu aboye the star There ig also an in . 

co-rotation radius, and the matter in the disk inside this region dication of Hne profile modulation b the ste llar rotation period 

free-falls onto the stellar photosphere following magnetic field , iBouvier et al Il2007h 
lines. When the plasma reaches the end of the accretion column, 
the impact with the stellar surface c reates a h ot spot, where its 

kineti c energy is thermalized (e.g., ICamenz ind 1990; iKoenigll However, the MA scenario is only partially successful in 

11991b . The gas is heated to temperatures of a 10 6 K in the shock, explaining the features observed in line profiles of CTTS. 

which then emits strongly in the X-rays. Most of the X-rays are Observational studies show evidence of outflows in these ob- 

then reabsorbed by the accretion column, a nd re-emitted as the jects, which are inferred by the blue-shifted absorption features 

blue and ultraviolet continuum excess (e.g.. lCalvet & Gullbririgl that can be seen in so me of the Balmer and Na D lines (e.g., 

|l998i) , which can be observed in CTTS. lAlencar & Basril UOQO). Bipolar outflows leaving these objects 

The magnetospheric accretion paradigm is currently strongly have been observed dow n to scales of =* 1 AU (e.gjTakami et al.1 

supported by observations. Magnetic field measurements in 5003; Appenzeller et al. 2005 ). On larg er scales, HST observa- 

tions of HH30 (e.g jBurrows et al.l 19961) could trace a bipolar jet 

* E-mail: styx@fisica.ufmg.br to within < 30 AU of the star. This suggests a greater complex- 
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ity of the circumstellar environment, in which the MA model is 
only a piece of the puzzle. 

The observed jets are believed to be highly collimated 
magneto-hydrodynamic disk winds which efficiently extract an- 
gular momentmn_and_gravitational energy from the accretion 
disk. iBlandford & Pavnd (Il982l hereafter BP) were the first 
to propose the use of a disk wind to explain the origin of 
jets from accreting disks aroun d a black hole, and soon after 
iPudritz & Norman ( 19831 fl986l) pro posed a similar mech anism 
for the origin of the protostellar jets. iBacciotti et al.l (120031) . us- 
ing high resolution spectro-imaging and adaptative optics meth- 
ods on protostellar jets, observed what seems to be rotational 
motion inside these jets, and also an onion-like velocity struc- 
ture where the highest speeds are closer to the outflow axis. 
Their work strongly corroborates the idea of the origin of jets 
as centrifugally driven MHD winds from extended regions of 
their accretion disks. A correlation between je ts, infrared excess, 
and the accretion pro cess has been observed (ICabrit et al.l l 1990: 
Hartigan et al.l 1 19951) . which indicates a dependence between 
the outflow and inflow of matter in CTTS. Disk wind theories 
show the existence of a scale relation between the disk accretion 
(M acc ) and the mass loss rates (M loss ) (Pelletier & Pudritz 1992). 
Observations have confirmed this relation (Hartmann 1998), and 
both have agreed that in CTTS the typical value for this rela- 
tion is Mi oss /M acc - 0.1. Another model widely used to explain 
the outflows and the protostellar j ets is called the X-wind model 
(IShuetal.l[l994t ICai et al.l 1200^) . which states that instead of 
originating in a wider region of the accretion disk, the outflow 
originates in a very constrained region around the so called 'X- 
point', which is at the Keplerian co-rotation radius of the stellar 
magnetosphere. 

Many developments in MHD time-de pendent numerica l 
models have been made in the last years. IZanni et all £2007) 
were able to reproduce the disk-wind lauching mechanism and 
jet formation from a magnetized accretion disk with a resistive 
MHD axisymmetric model showing the system evolution over 
tens of rotation periods. They were able to find a configuration 
in which a slowly evolving outflow leaving the inner part of the 
accretion ring (from ~ 0.1 to ~ 1 AU) was fo rmed, when using a 
high disk magnetic resistivity. Mu rphy et al.l d2010h showed that 
it was possible to form a steady disk wind followed by a self- 
contained super- fast-magnetosonic jet using a weakly magne- 
tized accretion disk. Their numerical solution remained steady 
for almost a thousand Keplerian orbits. Both m odels help to 
strengthen the disk wind scenario. Romanov a et al.l (12009). how- 
ever, showed that it is also possible to numerically simulate the 
launching of a thin conical wind, which is similar in some re- 
spects to an X-wind, considering only an axisymmetric stellar 
dipolar magnetic field around a slowly rotating star. A fast ax- 
ial jet component also appears in the case of a fast rotating star. 
Their conic al wind also shows som e degree of collimation, but, 
as stated by Roman ova et al.l d2009l) . this collimation may not be 
enough to explain the observed well-collimated jets. 

Radiative transfer models were used to calculate some of the 
line profiles that were observed in CTTS. The first of these mod- 
els that used the MA paradigm, hereafter the HHC model, was 
proposed by Hartmann, Hewett, & Calvet ( 1994), and was based 
on a simple axisymmetric dipolar geometry for the accretion 
flow, using a two-level atom approximation under the Sobolev 
"resonant co-moving surfaces" approximation (SA). The model 
was later improved by the addition of the full statistical equi- 
librium equations in the code ( Muzerolle. Calvet. & Hartmann! 
1998 ), and finally including an exact integr ation of the line pro- 
file (IMuzeroll e. Calvet. & Hartmann! 120 01). These models in- 



cluded only the magnetospheric and photospheric components 
of the radiation field and were partially successful in reproduc- 
ing the observed line strengths and morphologies for some of the 
Balmer lines and Na D lines. In these models, the stellar photo- 
sphere and the inner part of the accretion disk partially occult 
the inflowing gas, and these occultations create an asymmetric 
blue-ward or red-ward emission peak, depending on the angle 
between the system symmetry axis, and the observer's line of 
sight. The inflowing gas projected against hot spots on the stellar 
surfa ce produces t he IPC profiles that are sometimes observed. 

lAlencar et al.l (120051) demonstrated that the observed Ha, 
H/3, and Na D lines of RW Aur are better reproduced if a disk 
wind component arising from the inner rim of the accretion disk 
was added to these radiative transfer models, and it became clear 
that for more ac curate predictions, thes e models should include 
that component. iKurosawa et al. (2006, henceforth KHS), were 
the first to use a self-consistent model where all these compo- 
nents were included, and could reproduce the wide variety of the 
observed Ha line profiles. Still, their cold MHD disk wind com- 
ponent used straight field lines, and some of which violate the 
launching conditions required by BP, which state that for an out- 
flow to appear, its launching angle should be > 30° away from 
the rotational axis of the system. Those infringing field lines in 
the KHS model were always the innermost ones, which, as we 
will show here, are responsible for the bulk of the Ha line pro- 
files in the disk wind region. 

The BP launching condition assumes a cold MHD flow, 
which is not always the case. In a more general situation, the 
thermal pressure term becomes important near the base of the 
outflow, and should not be neglected during the initial accelera- 
tion process, and thus it is possible to have a stable steady solu- 
tion even if the launching angle exceeds 30°. The BP self-similar 
solution, although it is a very simple one, is able to describe 
mathematically and self-consistently the disk-wind launching 
mechanism and the further collimation of the disk wind to a 
jet by magnetic "hoop" stress. The self-similarity hypothesis 
becomes increasingly artificial farther away from the accretion 
disk, in the region where the jet is formed, but it remains valid 
near the accretion disk, where our calculations were performed. 

In this work, we have included a disk wind component in the 
HHC model, using a modified version of the BP formulation. We 
then studied the parameter space of the improved model, while 
trying to identify how each parameter affected the calculated Ha 
profiles. In Sects 2 and 3 we present the magnetospheric accre- 
tion, the disk wind and the radiative transfer models that were 
used as well as all the assumptions made to calculate our pro- 
files. The model results are shown in Sect. 4, followed by a dis- 
cussion in Sect. 5 and the conclusion in Sect. 6. A subsequent 
paper will focus on the modeling of the observed Ha lines for a 
set of classical T Tauri stars that exhibit different accretion and 
environmental characteristics. 



2. Model structure 

Our model structure comprises four components: the star, the 
magnetosphere, the accretion disk, and the disk wind. The cen- 
ter of the star is also the center of the model's coordinate sys- 
tem, and the star rotates around the z-axis, which is the system's 
axis of symmetry. We assume symmetry through the x-y plane. 
In our coordinate system, r is the radial distance from the cen- 
ter of the star, and 8 is the angle between r and the z-axis. The 
modeled region is divided into a grid system, where 1 /4 of the 
points are inside the magnetosphere, and the other 3/4 are part 
of the disk wind region. The star and accretion disk are assumed 
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Fig. 1. Not-to-scale sketch of the adopted axisymmetric geome- 
try shown in the poloidal plane. The gray half-circle denotes the 
stellar photosphere, the dark rectangle represents the disk, which 
is assumed to be opaque and to extend to the innermost magnetic 
field line. The dotted lines represent the magnetosphere and its 
purely dipolar field lines, and the dash-dotted lines are parabolic 
trajectories that represent the disk wind streamlines. The inner 
and outer magnetospheric field lines intersect the disk at r m \ and 
r mo , respectively, and the ring between rdi and r^ is the region 
where the disk wind is ejected. The inclination from the z-axis is 
given by 6, and §q is the disk-wind launching angle. 



to be infinitely thick in the optical region of the spectrum, and as 
such can be said to emit as blackbodies, and act as boundaries in 
our model. Due to symmetry arguments, it is only necessary to 
calculate the velocity, density, and temperature fields in the first 
quadrant of the model. These values are rotated around the z- 
axis, and are then reflected through the x-y plane. The emission 
from the accretion disk in CTTS is neglected, because its tem- 
perature is usually much lower than the photospheric tempera- 
ture of the star. A schematic diagram of the adopted geometry 
can be seen in Fig.Q] 

The stellar surface is divided into two regions: the photo- 
sphere and the accretion shock (the region where the accru- 
ing gas hits the stellar photosphere). Unless otherwise stated, 
the stellar parameters used in the model are those of a typical 
CTTS, i.e. radius (R*), mass (M») and effective temperature of 
the photosphere (T p h) are 2.0 R©, 0.5 M and 4000 K, respec- 
tively. When the gas hits the star, its kinetic energy is thermalized 
inside a radiating layer, and we consider that this layer emits as a 
single temperature blackbody. Because we assume axial symme- 
try, our accretion shock region can be regarded as two spherical 
rings near both stellar poles, and we as sume their effective tem - 
perature (T sh ) to be 8000 K, as used bv lMuzerolle et all J2001I) . 
Although T s h has been kept constant, it actually depends on the 
mass accretion rate. But it barely affects the Ha line profile and 
is only important when defining the continuum veiling levels. 



2. 1 . Magnetospheric model 

The basic assumption of this model is that accretion from the 
disk onto the CTTS is controlled by a dipole stellar magnetic 
field. This magnetic field is assumed to be strong enough to trun- 



cate the disk at some height above the stellar surface, and also to 
remain undisturbed by the ionized inflowing gas. The gas pres- 
sure inside the accretion funnel is assumed to be sufficiently low 
so that the infalling material free-falls onto the star. 

For the considered axisymmetric flow, the streamlines can be 
described by 



r = r m sin 6, 



(1) 



where r m corresponds to the point where the stream starts 
at the disk surfac e, or, in other words, where 6 - n/2 
dGhosh et al.l [1977). In this problem, we can use the ideal 
magneto-hydrodynamic (MHD) approximation, which assumes 
an infinite electrical conductivity inside the flow. According to 
this approximation, the flow is so strongly coupled to the mag- 
netic field lines that the velocity of the inflowing gas is parallel 
to the magnetic field lines, and thus the poloidal component of 
the velocity is 



3y l/2 (l- y y /2 v7 + (2-3y)z 



(4- 3V) 1 / 2 



(2) 



where y - r/r m - sin 0, m is the cylindrical radi al direction, and 
li p is th e gas free-fall speed, which according to lHartmann et al.l 
( 1994) can be written as 



2GM„ (R, R, 
R, 



1/2 



(3) 



This model constrains the infalling gas between two field 
lines, which intersect the disk at distances r mo , for the outermost 
line, and r m ;, for the innermost one (Fig.[T). The area filled by the 
accretion hot spots directly depends on these values, because it 
is bounded by the limiting angles 6>; and O , the angles where the 
lines crossing the disk at r m \ and r mo hit the star, respectively. We 
have set these values to be r m ;=2.2 R* and r mo = 3.0 R„, which are 
the same values used in the 'small /wide' case oflMuzerolle et al.l 
(2001). The corresponding accretion ring covers an area of 8% 
of the total stellar surface, which is on the same order of magni - 
tude as the observational estimates (Gullbring et al. 1998, 2000). 
Hartman n et al.l (11994) showed that the gas density (p) can be 
written as a function of the mass accretion rate M acc onto these 
hot spots as 



M„, 



-5/2 



P = 



4-3y 



47r(l/r ml -l/r mo )(2GM)'/ 2 V 1 -y 



(4) 



Martin (1996) presented a self-consistently solved solution 
for the thermal structure inside the dipolar magnetospheric ac- 
cretion funnel by solving the heat equations coupled to the rate 
equations for hydrogen. When his temperature structure was 
used by Muzerolle et al. (1998) to compute CTTS line profiles, 
their results did not agree with the observations. This is not the 
case with theor etical profiles bas e d on th e calculated temperature 
distribution by iHartmann et al.l (1 19941) . which we adopt here. 
The temperature structure used is computed assuming a volu- 
metric heating rate oc r~ 3 and bal ancing the ene rgy input with 
the radiative cooling rates used by IHartmann et a l. ( 1982). 

We used a non-rotating magnetosphere in this model. The 
rotational speed inside the magnetosphere is much lower com- 
pare d to the f ree-fall speeds inside the accretion columns. 
Muzerolle et al. (2001) showed that the differences between pro- 
files calculated with a rotating and a non-rotating magnetosphere 
are minimal if the stellar rotational speed is < 20 km/s. 
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2.2. Disk wind model 

We adopted a modified BP model for the disk wind. This model 
assumes an axisymmetric magnetic field treading the disk, uses 
the ideal MHD equations to solve the outflowing flux leaving 
the disk, and assumes a cold flow, where the thermal pressure 
terms can be neglected. In the BP model, the wind solution is 
self-similar and covers the whole disk; in this work, instead, we 
use a self-similar solution only inside the region limited by the 
disk wind inner and outer radii. This solution can be specified by 
the following equations: 



V(pv) = 0, 



p(v • V)u z 



~ P Jz 



-GM t 



(m 2 +z 2 y 2 

kB 

v = - — + (w x r), 
Anp 

1 2 GM t 

6 ~ 2 V ~ (m 2 +z 2 y/ 2 



(5) 

1 r)R^ 1 

8^ + ^ (B ' V) ^ (6) 



wfii 



/ 



■utBa, 



VJV$ 



(7) 
(8) 
(9) 



where e, the specific energy, /, the specific angular momentum, 
k, the ratio of mass flux to magnetic flux, and to, the angular 
velocity of the magnetic field line, are constants of motion along 
a flow line (Mestel 1961, BP). The variables r, v, B are the radial, 
the flow velocity, and the magnetic field vectors in cylindrical 
coordinates respectively. If we solve the problem for the first 
flux line, the solution can be scaled for the other lines with the 
scaling factors introduced in BP: 



(m^(x), 0, mox) , 



(10) 



GM t 

Co 
GM t 



mi 
P = PoVix) 



1/2 









<-/ r 3/2 
= Pl \~) ^> 



B = Bo (M*).W. **(*)) 

/ \-5/4 

X = Z/tDq, 



(13) 
(14) 



where a zero subscript indicates evaluation at the disk surface 
(X-Q, £=1), a subscript 1 means that the quantity is evalu- 
ated at the fiducial radius m\, and a prime denotes differenti- 
ation with respect to x- The functions %(x), /Of). g(x)> 7 l(x)^ 
b w (x) , bsix), b z (y) represent th e self-similar solution of the 
wind. Blandfo rd & Pavna (1 19821) introduced the dimensionless 
parameters 



(15) 
(16) 
(17) 



(GMJmoY 
I 



(GM,57 )" 2 

r/2 xi/2 (GMJmo) 1 ' 2 



k(l + {' 2 ) 



Bn 



with e, A, k and & the constant parameters that define the so- 
lution. With the help of Eqs. (fT0li-(fT7l). and knowing that the 



Keplerian velocity is a> — (GM /s^) 1 ^, one can find after some 
algebraic manipulation (see BP lSafieril 993) the following quar- 
tic equation for /Of), 



T 

where 



f 2 U 



(A - f)m 



&\-m) 



47T P (u 2 m + VJ) 

m - ; n — > 

Bl + Bl 

= KSfJ 

is the square of the (poloidal) Alfven Mach number, and 

T = f + 2S - 3, 

s = (t 2 +x 2 r l/ \ 
u = (i+r 2 ). 
J = t-xt- 



(18) 

(19) 
(20) 



(21) 
(22) 
(23) 
(24) 



The flow trajectory solution is a self-similar solution that 
always crosses the disk at £ = 1 . If this self-similar solution, 
which is given by ;t(£), is expanded into a polynomial function 
around the disk-crossing point using Taylor's theorem, and we 
consider only the first three terms of the expansion, we will have 
a parabolic trajectory describing the outflow. This parabolic tra- 
jectory is given by 



X(0 = a s +H + c, with c - -(a + b), 



(25) 



This solution respects the launching condition, which states that 
the flow-launching angle #o (see Fig.[TJ should be less than -& c = 
60° for the material to accelerate magneto-centrifugally and be 
ejected as a wind from the disk. Equation (TT~8T > can be rewritten 
as a quartic equation for m, 



fVm 4 - 2fVm 3 - [f(T -V)-(A- ff\ m 2 + 



(11) 

(12) where V = 



+2fTm + fT 



U 



(k?J) 2 



(26) 



(27) 



The solution of Eq. ( 1261 ) gives four complex values of m for 
each pair of points (x, £) inside the trajectory. The flow must 
reach super- Alfvenic speeds, for a strong self-generated toroidal 
field to appear, collimatin g the flux, and thus producing a jet 
dBlandford & Pavnall982h . Of the four solutions of this quartic 
equation, there are only two real solutions that converge on one 
another at the Alfven point, a sub-Alfvenic and a super- Alfvenic 
one. The sub-Alfvenic solution reaches the Alfven speed and 
then the solution starts decreasing again, thus remaining a sub- 
Alfvenic solution before and after the Alfven point. The same 
happens with the super- Alfvenic solution, which remains super- 
Alfvenic before and after the Alfven point. The chosen solution 
must be a continuous real solution, and must also increase mono- 
tonically with x, because the wind is constantly accelerated by 
the magnetic field. This solution is obtained by using the conver- 
gent sub-Alfvenic solution before the Alfven radius, and then the 
super- Alfvenic one. 

The launching condition imposes that & > cot# c , and it 
helps constrain the set of values of a, b, and c that can be used in 
equation ( fZSb . Then we can write 



fo = 



1 



2a + b 



> COt ft c 



(28) 
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If the values of parameters k and A are also given, the flow can 
be fully described by the following equations: 



m 

i ~ IdTr 

/ p = /(1+£' 2 ) I/2 , 
% 2 -mA 

8 ~ £(l-m)' 
(€fJ)o _ wo 



(29) 
(30) 
(31) 

(32) 



where 77 and / p are the self-similar dimensionless density and 
poloidal velocity, respectively. Knowing the gas density p, the 
poloidal velocity v p , and the scaling laws, it can be shown that 
the mass loss rate Mi oss is given by 



M loss = 27rp /po \GM t ml) In — , 



(33) 



where r^ and rd are the inner and outer radii in the accreting 
disk that delimits the wind ejecting region (see Fig.Q]). With this 
equation it is possible to calculate the value of po, which is nec- 
essary to produce a Mi oss , when the wind leaves the disk between 
r di and r do . 

The temperature structure we used inside the disk wind is 
similar to the one used in the magnetospheric accretion funnel. 
There are some studies about the tem perature structure inside 
this region dSafierl 1 993b [Panodou et al. 2010), and we are going 
to discuss them in Sect. 6. 

The parameters used to calculate the self-similar wind solu- 
tion were initially: a=0.43, b=-0.20, .1=30.0, k=0.03. These val- 
ues of a and b were chosen to produce an intermediate disk- wind 
lauching angle #o-33°, while for A and 6 these values are the 
same as the ones used by BP in their standard solution. We later 
changed some of these values to see how they influence our cal- 
culated Ha line profiles. The disk wind inner radius rdi=3.01 R t 
is just slightly larger than r mo , and we changed the value of the 
disk wind outer radius rd from 5 R t to 30 R* to see how this 
affected the calculated profiles. To calculate po we used a base 
value of Mioss = 1 0~ 9 M yr" ' =0. 1 M acc . 

3. Radiative transfer model 

In order to calculate the Ha line, it is necessary to know the 
population of at least levels 1 to 3 of the hydrogen atom. The 
model used for the hydrogen has three levels plus continuum. 
The high occupation number of the ground level, owing to the 
low temperatures in the studied region, leads to an optically thick 
intervening medium to the Lyman lines. Using this, it can be as- 
sumed that the fundamental level is in LTE. The other two levels 
plus ionization fraction were calculated with a two-level atom 
approximation (see Mihalas 1978, chapter 11). The fundamental 
level population is used to better constrain the total number of 
atoms that are in levels 2 and 3, and is necessary to calculate the 
opacity due to H ion. 

The line source function is calculated with the Sobolev ap- 
proximation method (Rvbicki & Hummer 1978; Hartman n et alj 
1994). This method can be applied when the Doppler broaden- 
ing is much larger than the thermal broadening inside the flow. 
Thus, an atom emitting a specific line transition only interacts 
with other atoms having a limited range of relative velocities, 
which occupy a small volume of the atmosphere. In the limit of 
high velocity gradients inside the flow, these regions with atoms 
that can interact with a given point additionally become very 



thin, and it can be assumed that the physical properties across 
this thin and small region remain constant. These small and thin 
regions are called "resonant surfaces". Each point inside the at- 
mosphere then only interacts with a few resonant surfaces, which 
simplifies the radiative transfer calculations a lot. The radiation 
field at each point can be described by the mean intensity 



/ V = [l- J 8(r)]5(r)+ J 8 c (r)7 c +F(r), 



(34) 



where ft and /3 C are the scape probabilities of local and contin- 
uum (stellar photosphere plus accretion shock, described by I c ) 
radiation, S is the local source function, and F is the non-local 
term, which takes into account only the contributions from the 
resonant surfaces. A more detailed description of each of these 
terms has been given in Hartm ann et al.1 (1994). 

For a particular line transition, the source function is given 
by 



S u i = 



2hv 3 



Nigu 

N u gi 



- 1 



(35) 



where Ni, N u ,gi and g u are the population and statistical weights 

of upper and lower levels of the transition respec tively. 

Th e flux is calculated as described by iMuzerolle et al.l 
(1200 lh . where a ray-by-ray method was used. The disk coor- 
dinate system is rotated to a new coordinate system, given by 
the cartesian coordinates (P, Q, Z), in which the Z-axis coincides 
with the line of sight, and each ray is parallel to this axis. The 
specific intensity of each ray is given by 



/-T(Zo) 
-Tot + I 

Jt{— do) 



/ v = / oe - T '°>+ I S v {r')Q- T &Y 
•(-00) 



(36) 



where Iq is the incident intensity from the stellar photosphere or 
the accretion shock, when the rays hit the star, and is zero oth- 
erwise. S V (Z) is the source function at a point outside the stellar 
photosphere at frequency v, and r tot is the total optical depth 
from the initial point, Zo, to the observer (-00), and r(-oo) = 0. 
Zo can be either the stellar surface, the opaque disk, or 00. The 
optical depth from a point at Z to the observer is calculated as 



T(Z) 



U—Oi 



[ Xc (Z')+Xi(Z')]dZ', 



(37) 



in which Xc and y\ are the continuum and line opacities. The 
continuum opacity and emissivity contain the contributions from 
hydrogen bou nd-free and fre e-free emission, H", and electron 
scattering (see lMihalaslll978l) . 

This model uses the Voigt profile as the profile of the emis- 
sion or absorption lines that were modeled. Assuming a damp- 
ing constant F, which controls the broadening mechanisms, 
a - F/47tAvd, v = (v - vo)/Avd, and Y - Av/Avq, where Avd is 
the Doppler broadening due to thermal effects, and vo the line 
center frequency, the Voigt function can be written as 



a r e- y2 

H(a,v)=- -= -dY. 

x J-00 (v - Y) 1 + a 1 



(38) 



The line opacity is a function of the Voigt profile, and can be 
calculated by 



x l/2 ql I gjUi \ 

X\ = fijTij 1 H(a, v), 

m e c ■ \ gill; ' 



(39) 



where fy, n,, ttj, gi, and gj are the oscillator strength between 
level i and j, the populations of the jth and j'th levels and the 
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Table 1. Model parameters used in the default case. 



Stellar Parameters 


Values 


Stellar Radius 


(R.) 


2.0 R 


Stellar Mass 


(M.) 


0.5 M Q 


Photospheric Temperature 


(T ph ) 


4000 K 


Magnetospheric Parameters 


Values 


Accretion Shock Temperature 


(T sh ) 


8000 K 


Mass Accretion Rate 


(M acc ) 


lO^Moyr" 1 


Magnetosphere Inner Radius 


Omi) 


2.2 R» 


Magnetosphere Outer Radius 


(/mo) 


3.0 R» 


Disk Wind Parameters 


Values 


Inner Radius 


OdO 


3.01 R, 


Outer Radius 


Odo) 


30.0 R« 


Mass Loss Rate 


(M loss ) 


0.1M acc 


Launching angle 


(fto) 


33.4° 


Dimensionless Angular Momentum 


w 


30.0 


Mass Flux to Magnetic Flux Ratio" 


(*) 


0.03 



"- Dimensionless parameter 

degeneracy of the /th and /'th levels, respectively. The electron 
mass and charge are represented by m e and q e . The damping 
constant T contains the half-width terms for radiative, van der 
Waals, and Stark br oadenings, and has been parameterized by 
IVernazza etail dl973h as 



T - C ra H + C 



\10 16 cm- 3 /\5000K/ 



(40) 



+C Sl 



tark 



/ N e \ 2/3 

\10 12 cm- 3 / ' 



where C ra d, C v dw and Cstark are the radiative, van der Waals, and 
Stark half-widths, respectively, in A, Nhi is the number density 
of neutral hydrogen, and N e the number density of electrons. 
Radiative and Stark broadening mechanisms are the most impor- 
tant for the Ha line in our modeled environment. Stark broaden- 
ing is significant (or dominant) inside the magnetosphere, due 
to the high electron densities in this region. However, inside the 
disk wind, N e becomes very low, and Stark broadening becomes 
very weak compared to radiative broadening. These effects are 
small compared to the thermal broadening inside the studied re- 
gion. 

The total emerging power in the line of sight for each fre- 
quency (P Y ) is calculated by integrating the intensity I v of each 
ray over the projected area dPdQ. The emerging power in the 
line of sight depends only on the direction. Thus, instead of us- 
ing flux in our plots, which also depends on the distance from 
the source, we present our profiles as power plots. 

4. Results 

As a starting point, we used the parameters summarized in 
Table Q] to calculate Ha line profiles. We used a default in- 
clination /=55° in these calculations. For Mi oss , r^ and r^ 
as given by Table [TJ we calculated with Eq. (l33l a fiducial 
density pa-A.l x lCT 12 gcirr 3 . The set of maximum tempera- 
tures inside the magnetosphere and the disk wind region em- 
ployed in the model ran ges from 6000 K to 10 000 K, as used by 
iMuzerolle et al.l d200ll) . As an example, Fig. [2] shows the tem- 
perature profiles for the case when the maximum temperature 
inside the magnetosphere, r mag MAx - 8000 K, and the max- 
imum temperature inside the disk wind, r w j„d,MAx — 9000 K. 
Figure|2]shows some of the streamlines to illustrate the tempera- 
ture behavior along the flow lines. Inside the magnetosphere, we 



can see that there is a temperature maximum midway between 
the disk and the star, and the maximum temperature is higher for 
streamlines starting farther away from the star. Inside the disk 
wind it is the opposite, the maximum temperature along a flux 
line starts higher near the star, and decays with receding launch- 
ing point from the star. Also for the nearest flux lines the temper- 
ature reaches a maximum and then starts falling inside the disk 
wind, and for more extended launching regions the temperature 
seems to reach a plateau and then remains constant along the 
line. The right plots in Fig.[2]show the density profiles along the 
streamlines inside the magnetosphere and disk wind. The tem- 
perature law used assumes that T oc N H 2 r^, where Nh is the to- 
tal number density of hydrogen, which explains the similarities 
between the temperature and density profiles. Each streamline 
is composed of 40 grid points. These lines reach a maximum 
height of 32 R* in the disk wind. 

Figure [2] also shows the poloidal velocity profiles inside the 
magnetosphere and the disk wind. Inside the magnetosphere, the 
velocity reaches its maximum when the flow hits the stellar sur- 
face, and it is faster for lines that start farther away from the star. 
Inside the disk wind, the outflow speed increases continuously 
and monotonically, as required by the BP solution, and this ac- 
celeration is fastest for the innermost streamlines. We can see 
that the terminal speeds of each streamline inside the magne- 
tosphere have a very low dispersion and reach values between 
^230 kms 1 and ^260 kms" 1 . Inside the disk wind, this veloc- 
ity dispersion is much higher, with the flow inside the innermost 
streamline reaching speeds faster than 250 km s _1 , while the flow 
inside the outermost line shown in Fig. [2]- which is not identi- 
cal with the outermost line we have used in our default case - 
is reaching maximum speeds around ^ 70 km s _1 . Note that even 
considering a thin disk, our disk wind and magnetosphere so- 
lutions do not start at z=0, and the first point in the disk wind 
streamlines is always at some height above the accretion disk, 
which is the reason why the poloidal velocities are not starting 
from zero. 

The Ha line profiles have a clear dependence on the tem- 
perature structures of the magnetosphere and the disk wind re- 
gion. Figure [3] shows plots for different values of r mag; MAX, and 
each plot shows profiles for different values of r w j n d,MAx, which 
are shown as different lines. The black solid lines are the pro- 
files with only the magnetosphere's contribution. One can see 
that in the default case, most of the Ha flux comes from the 
magnetospheric region, which is highly dependent on the mag- 
netospheric temperature. In all plots there is no noticeable differ- 
ence between T^^max - 7000 K (red dot-dot-dot dashed line) 
and r w ; n d,MAX - 8000 K (green dash-dotted line), and the disk 
wind contribution in these cases is very small, only appearing 
around the line center as an added contribution to the total flux, 
while there are no visible contributions in the line wings. When 
T'wind.MAx is around 9000 K, a small blue absorption component 
starts appearing, which becomes stronger as both the magneto- 
spheric and disk wind temperatures increase. The same happens 
for the added disk wind contribution to the flux around the line 
center. 

Figure [4] shows the variation of the Ha line with the mass 
accretion rate (M acc ). In all plots, the mass loss rate was kept 
at 10% of the mass accretion rate, and the maximum temper- 
ature inside the magnetospheric accretion funnel was around 
9000 K. For M acc =lO _9 M yr- 1 , the plot shows that the line 
does not change much when the disk wind temperature varies 
from 6000 K to 10 000 K, because all lines that consider the disk 
wind contribution are overlapping. There is an overall small disk 
wind contribution to the flux that is stronger around the line 
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Fig. 2. Temperature (left), poloidal velocity (middle), and density (right) profiles along streamlines inside the magnetosphere (top) 
and along streamlines inside the disk wind region (bottom). In all plots the x-axis represents the system cylindrical radius in 
units of stellar radii R». The maximum temperature inside the magnetosphere is T^^max — 8000 K, and inside the disk wind it 
is ^wind.MAX - 9000 K. Only some of the streamlines are shown in the plots, and they are the same ones in the left, middle, and right 
plots. 



center. In this case, the disk wind contribution to the line cen- 
ter increases slightly as the temperature becomes hotter than 
10 000 K. If the temperature becomes much higher than that, 
most of the hydrogen will become ionized and no significant ab- 
sorption in Ha is to be expected. When the mass accretion rate 
rises to M acc =lO _8 M yr _1 , the Ha flux becomes stronger, and 
the line broader. There is no blue-ward absorption feature until 
the temperature inside the disk wind reaches about 9000 K (see 
Fig. [3]), after that the blue-shifted absorption feature becomes 
deeper as the temperature rises. Finally, for a system where 
M acc =10 _7 M Q yr _1 , the line becomes much broader, and even 
for maximum temperatures inside the disk wind below 8000 K, 
a blue-shifted absorption component can be seen. This absorp- 
tion becomes deeper with increasing disk wind temperature. In 
this last case, the disk wind contribution to the flux is on the same 
order of magnitude as the contribution from the magnetosphere, 
which does not happen for lower accretion rates. 

When po is kept constant and the outer wind radius varies, 
then Mioss also changes. When that is done, it is also possible 
to discover the region inside the disk wind which contributes 
the most to the line flux. Figure [5] shows how the line changes 
with rdo, while po remains constant. Where M acc =10 _8 M Q yr _1 
(Fig. |5h), we can see that the line profiles for rd o =30.0R» and 
7"do=10.0R* are almost at exactly the same place, which means 
that in this case most of the contribution of the disk wind comes 
from the wind that is launched at r < 10 R„. And as r^ becomes 
even smaller, the blue-shifted absorption feature begins to be- 
come weaker and the line peak starts to decrease, until the only 



contributing factor to the line profile is the magnetospheric one. 
For M acc =10~ 7 M yr~' (Fig.[5t>), a similar behavior can be seen, 
but in this case the plots show the loss of a substantial amount 
of flux when rdo goes from rd o =30.0R» to rd o =10.0R», an ef- 
fect mostly seen around the line peak. The same happens at the 
blue-shifted absorption feature, which begins to weaken when 
/"do ^15.0 R„, and weakens faster as r do decreases. In both cases, 
most of the Ha flux contribution from the wind comes from the 
region where the densities and temperatures are the highest, and 
this region extends farther away from the star for higher mass 
loss rates, as expected. But even for high values of Mi oss , all ab- 
sorption comes from the gas that is launched from the disk in a 
radius of at most some tens of R, from the star. In both cases, it 
seems that a substantial part of the disk wind contribution to the 
Ha flux around the line center comes from a region very near the 
inner edge of the wind, because we still can see a sizeable disk 
wind contribution around the line center even when rd =3.5 R*. 

If instead Mi 0S8 is kept constant while both r do andpo change, 
it is possible to infer how the size of the disk wind affects 
the Ha line as shown in Fig. [6] In the default case (Fig. [6^, 
M^lO^Moyr 1 ), with r mag ,MAx « 9000 K, and r wind , M AX - 
10 000 K, and changing ^ from 30 R* to 5 R* while also chang- 
ing the densities to keep the mass loss rate constant, there are 
almost no variations around the line peak, as all the profiles 
with the disk wind contribution overlap in this region. Some 
differences can be seen at the blue-shifted absorption compo- 
nent, which becomes deeper as the disk wind size decreases 
and the density increases. For a higher mass loss rate (Fig. [6J}, 
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Fig. 3. Ha profiles for the default case with values of r mag MAX ranging from 7500 K to 10000K and values of r w j„d,MAx ranging 
from 7000 K to 10 000K. The maximum magnetospheric temperature r mag ,MAX in eacn plot is the same, and each different line 
represents a different maximum disk wind temperature T^^max'- 7000 K (black solid line), 8000 K (grey solid line), 9000 K (black 
dotted line) and 10 000 K (black dashed line). The dash-dotted dark grey lines in each plot show the profiles calculated with only 
the magnetospheric component (no disk wind). The black and grey solid lines completely overlap each other. 



M acc =10~ 7 M Q yr _1 ), the line peak decreases when r& goes from 
30 R* to 20 R„, while the absorption feature barely changes. But 
as rjo diminishes even further, the line flux as a whole begins 
to become much stronger again, while the absorption appar- 
ently becomes lower. Thus, it seems for moderate values of M acc 
that the size of the disk-wind launching region does not have 
a great influence on the Ha line flux, with only some small 
differences in the depth and width of the line profile absorp- 
tion feature. However, for higher M acc , the profiles start to vary 
rapidly when rdo < 20.0 R», and the line flux increases faster as 
rdo decreases. In spite of these differences, we noticed that in 
both cases this increase in the Ha flux only becomes important 



when po > 10 10 gcm 3 , 



which happens for rdo 5 7.5 R, when 



M acc =10 M yr , and for r do <4R», which is out of the plot 
range in Fig. [6^, when M acc =10 _8 M Q yr _1 . This indicates that 
the densities inside the disk wind region are very important to 
define the overall Ha line shape, if po > 10~ 10 gcrrT 3 . Still, if 
po is smaller, its effect, if any, is smaller and is mostly present 
around the absorption feature. 

So far we used A=30 and k-0.03 in all above mentioned 
cases. The A parameter is the dimensionless angular momentum, 
and together with e, the dimensionless specific energy, is con- 
stant along a streamline. Equations ([8]), (O, d!5t . and ( fT6b show 
that a variation in A will produce a change in the velocities along 
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Fig. 4. Ha profiles for different values of M acc : (a )1(T 9 M yr~' , 
(b)lO~ 8 M yr~' and (c)10~ 7 Moyr _1 . The maximum magneto- 
spheric temperature used is r mag MAx-9000K, and the values 
used for r w i n( j,MAx are 6000 K (black thin solid line), 8000 K 
(dark grey thick dashed line) and 10000K (light grey thin 
dashed line). The light grey thin solid lines are the profiles 
where only the magnetosphere is considered. In all cases, the 
Mioss~0.1M acc . 



the streamline. This change in A will also produce a different 
disk wind solution, and consequently it is necessary to change 
the value of po to keep the mass loss rate the same. The effect of 
A in the line profile can be seen in Fig. |7h for the default model 
parameters and a constant mass loss rate. It shows that as the 
value of A is increased, the absorption feature becomes deeper 



and closer to the line center, which means a slower disk wind, 
and also an increase in the line flux around its center. Figure |7j5 
illustrates this behavior and shows how the poloidal velocity pro- 
file changes when the value of A is varied. 

The Ha flux varies in a similar manner when the value of k, 
the dimensionless mass flux to magnetic flux ratio, is changed, 
which is illustrated by Fig. [7J;. As k increases, the mass flow 
becomes slower (see Eqs.[TTIand|29j>, and the absorption feature 
also comes closer to the line center. The higher the value of k, 
the deeper is the blue-shifted absorption and the stronger is the 
flux at the line center, similar to the case when A varies. There is 
also a large variation of the poloidal speed inside the flux when 
k changes, which is shown by Fig. [TJl. Both Figs. [7J3 and [Vji 
confirm the behavior seen in Figs.|7^ and[7J;. An increase in the 
values of A or k decreases the poloidal velocities and the velocity 
dispersion inside the disk wind, causing the absorption feature 
to approach the line center. The lower velocity dispersion causes 
the line opacity to be raised around the absorption feature, which 
strengthens it. This increase in A or k also requires an increase in 
the densities inside the outflow, if the mass loss rate is to be kept 
constant, and this effect also enhances the absorption feature and 
the line flux around its center. 

The disk-wind launching angle §0 also produces a variation 
in the Ha line profile. To change -&q, it is necessary to change 
the coefficients a and b of Eq. ( f25l ). Here, b is kept constant, and 
only the coefficient a is varied. A change in the coefficients a and 
b produces another disk wind solution, and then it is necessary 
to recalculate po to hold the mass loss rate constant. In Fig. [8] 
■&o varies from 14.6° to 55.6°, which is almost at the launching 
condition limit, and the value of a used in each case is indi- 
cated in the caption. Figure [8] shows that the absorption feature 
moves closer to the line center as the launching angle steepens. 
Again, as with the cases when we changed k and A, a variation in 
the launching angle produces a variation in the launching speed, 
and this variation is shown at the bottom panel of Fig. [8] The 
launching speed is slower for higher launching angles, and this 
increases the densities inside the wind. When #o=14.6°, the case 
with the fastest launching speed, there is no noticeable differ- 
ence between the profile with the disk wind contribution and the 
one with only the magnetosphere contribution. This is also the 
case with the lowest disk wind densities, which are so low that 
there is no visible absorption, and just a small disk wind con- 
tribution around the center of the line. The case with i?o=33.4° 
is the default case used in all the prior plots. When the launch- 
ing angle becomes steeper, the disk wind becomes denser, and 
then we have a larger disk wind contribution to the profile and a 
deeper blue-shifted absorption feature. Again, the lower veloc- 
ity dispersion inside the disk wind also helps to enhance the disk 
wind contribution to the Ha when the launching angle becomes 
steeper. 

One of the most important factors defining the line profiles 
is the system inclination from the line of sight, i. Figure [9] shows 
how the Ha line changes with i, when M acc =10 _8 M o yr _1 , 
7mag,MAx=8000 K and r w i n d,MAx=10 000 K. It is obvious that the 
profile intensity is stronger for lower inclinations, and this has a 
series of reasons. First, because the solid angle containing the 
studied region is larger for lower values of i, and this leads to a 
higher integrated line flux over the solid angle. The opacity in- 
side the accretion columns is lower if we are looking at them 
at a lower inclination angle, and then we can see deeper in- 
side, which increases the line flux. Finally, the shadowing by the 
opaque accretion disk and by the star makes the lower accretion 
columns less visible with a higher i, which also decreases the to- 
tal flux. The Ha photon that leaves the magnetosphere or the star 
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Fig. 5. Ha profiles for (a) M acc =lCT s M yr~ , and 
(b) M acc =lO- 7 M yr- 1 , with 7 mag>MA x=9000K and 
7wind,MAx=10000K. In all plots the value of po is kept constant 
at a value which makes Mi oss x:0.1 M acc , when rd o =30.0R». The 
black thin solid lines are the line profiles with only the magneto- 
sphere; the thick darkgrey solid lines are the line profiles when 
magnetosphere and disk wind are considered and rd o =30.0R*. 
The other curves show the line profiles when rd o =10.0R» 
(thin black dashed lines), rd o =5.0R* (thin black dotted line), 
rdo = 4.0 R* (dark grey thin dashed line) and r& - 3.5 R* (light 
grey thin solid line). 



crosses an increasingly larger portion of the outflow as /becomes 
steeper, which strengthens the blue-shifted absorption feature. 
Another effect that happens when the inclination changes is that 
the disk wind velocities projected on the line of sight decrease 
as the inclination becomes steeper, and make the blue-shifted 
absorption feature approach the line center, as is illustrated by 
Fig-El The projected velocities increase until i reaches a certain 
angle, which depends on the disk wind geometry, and then as the 
inclination keeps increasing, the projected velocities start to de- 
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Fig. 6. Ha profiles as in Fig. [5] but now the density po is 
varied to keep the mass loss rate at 0.1 M acc , while the disk 
wind outer radius r^o changes from 30.0 R* to 5.0 R*, and 
each line represents a different rdo and a different po. In (a), 
r do =30.0R* and p =4.71xl0~ 12 gcnr 3 (dark grey thick solid 
line), r do =20.0R* andp =5.72xlO _12 gcrrT 3 (black thin dashed 
line), r do =10.0R* and p =9.02xl0 "" 
dotted line), rdo 

thin dashed line), r do =5.0R« and p =2.13xl0" n gcm 3 (light 
grey thin solid line). In (b), the lines have the same value of rdo, 
but their po are 10 times higher than in (a). The black thin line 
represents the line profile with only the magnetospheric compo- 
nent. 



crease again. The disk wind contribution to the red part of the Ha 
profile becomes weaker as i decreases, because there is barely a 
contribution to the red part of the profile when z'=15°, while this 
red-shifted contribution seems to be much stronger when /=60° 
and /=75°. If the disk-wind launching angle, §o, is smaller than 
i, part of the disk wind flux will be red-shifted. This red-shifted 
region inside the outflow becomes larger as i increases, and this 
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increases the red-ward disk wind contribution to the total line 
flux. 

The system inclination also affects the continuum level in 
the line profile. At lower inclinations, the upper accretion ring, 
with its much higher temperature than the stellar photosphere 
is much more visible than at higher inclinations. Moreover, at 
higher inclinations most of the stellar photosphere is covered by 
the accretion columns. Then, in a system with low inclination, 
the continuum level is expected to be significantly stronger than 
in systems with a high inclination. Analyzing Figs. [4] and [9] to- 
gether, we can see that this variation in the continuum is really 
significant. Figure [4^ shows exactly where the continuum level 
is when j=55°, and comparing it with Fig. [9] it is possible to note 
an increase between 30% and 50% in the continuum level when 
the inclination changes from 75° to 15°. 

Many of the line profiles calculated by our model have some 
irregular features at their cores. These features become more ir- 
regular as the mass accretion rate and temperature increase. A 
very high line optical depth inside the magnetosphere due to the 
increased density can create some self-absorption in the accre- 
tion funnel, resulting in these features. Also, the Sobolev ap- 
proximation breakes down for projected velocities near the rest 
speed. Under SA, the optical depth between two points is in- 
versely proportional to the projected velocity gradient between 
these points. Thus, in subsonic regions like in the base of the 
disk wind, or in regions where the velocity gradient is low, SA 
predicts a opacity that is much higher than the actual one, which 
creates the irregularities seen near the line center that are present 
in some of the profiles. If instead we do not assume SA, our 
profiles would be slightly more intense and smoother around the 
line center. However, the calculation time would increase con- 
siderably. 

5. Discussion 

The analysis of our results indicates a clear dependence of the 
Ha line on the densities and temperatures inside the disk wind 
region. The bulk of the flux comes mostly from the magneto- 
spheric component for standard parameter values, but the disk 
wind component becomes more important as the mass accretion 
rate rises, and as the disk wind temperature and densities become 
higher, too. For very low mass loss rates (Mi oss = 10~ 10 M Q yr -1 ) 
and disk wind temperatures below 10000K, there is a slight 
disk wind component that helps the flux at the center of the 
line only minimally. No absorption component can be seen even 
for temperatures in the wind as high as 17 000K, which is ex- 
pected because in such a hot environment most of the hydro- 
gen will be ionized. Thus, only the cooler parts of the disk wind 
where the gas is more rarified would be able to produce a min- 
imal contribution to the He flux. The fiducial density in this 
case is po = 4.7 x 10~ 13 gcirT 3 , when the disk- wind launching 
region covers a ring between 3. OR* (~ 0.03 AU) and 30. OR* 
(« 0.3 AU). Some studies suggest that the outer launching radius 
for the atomic component of the disk wind i s between 0.2 — 3 AU 
(lAnderson et al.ll2003b IPesenti et al.ll2004l iFerreira et al.ll2006t 
iRav et a02007h . which is much larger that the values of the disk 
wind outer radius used here. This means that the densities in- 
side the disk wind should be even lower than the values we used, 
but as the fiducial density falls logarithmically with r& Q (Eq.[33l). 
the density would not be much lower. The density for low mass 
loss rate (10~ 10 M yr -1 ) is already so low that a further decrease 
in its value would not produce a significant change in the disk 
wind contribution to the profile. However, the decreased den- 
sity would create an environment where a blue-shifted absorp- 



tion in the Ha line would be much more difficult to be formed. 
So, for M acc < 10~ 9 M G yr -1 , the disk wind contribution to the 
total Ha flux is negligible, w hich is consistent with observations 
(e.g. iMuzerolle et alJtzOOl 12005). A blue-shifted absorption in 
this case is meant to appear if the disk wind outer radius is much 
closer to the star than the observational studies suggest. 

For cases with higher values of M acc , and thus higher den- 
sities inside the disk wind region, it is possible to see a transi- 
tion temperature, below which no blue-shifted absorption fea- 
ture can be observed (Fig. @J, This temperature is between 
8000 K and 9000 K when M acc = 10~ 8 M yr" 1 , between 6000 K 
and 7000 K when M acc = lO _7 M yr _I , and it might be even 
lower than 6000 K for cases with very high accretion rates 
(M acc > 10 _7 M o yr _1 ). This transition temperature will directly 
depend on the mass densities inside the disk wind. Figure|4]also 
shows that the overall added disk wind contribution to the flux 
becomes more important as the accretion rate becomes stronger, 
due to a denser disk wind. It is then possible to imagine some 
conditions that would make the disk wind contribution to the 
Ha flux even surpass the magnetospheric contribution: a very 
high mass accretion rate (M acc > 10 _7 M o yr _1 ), or a very large 
region inside the disk wind with temperatures ~ 9000 K, or 
rdo <k 30 R«, or even a combination of all these conditions. 

The density inside the outflow alone is a very important fac- 
tor in the determination of the disk wind contribution to the Ha 
total flux. The blue-shifted absorption feature that can be seen 
in many of the profiles is directly affected by the disk wind den- 
sities, which change their depth and width as po changes. Yet, 
the other parts of the line profile seem barely to be affected by 
the same variation in po, if po < 10~ 10 gcrrr 3 , and the mass 
loss rate and the disk wind solution remain the same. But when 
Po <: 10~ 10 gcirr 3 , the variation in the line profile begins to 
adopt another behavior, and the line as a whole becomes depen- 
dent on the disk wind densities, which becomes stronger as po 
increases even more. Densities higher than that transition value 
would only occur, again, if the mass loss rates are very high, or if 
rdo is just a few stellar radii above the stellar surface in the case 
of moderate to high mass accretion rates. 

Another interesting result that we obtained shows that even 
if the disk wind outer radius ra is tens of stellar radii above 
the stellar surface, most of the disk wind contribution to the Ha 
flux around line center comes from a region very near the disk 
wind inner edge. This happens because of the higher densities 
and temperatures in the inner disk wind region if compared to 
its outer regions, which makes the inner disk wind contribution 
much more important to the overall Ha flux than the other re- 
gions. This result suggests that the more energetic Balmer lines 
should be formed in a region narrower than the one where the Ha 
line is formed, but a region that also starts at r^. However, the 
blue-shifted absorption feature seems to be formed in a region 
that goes from the inner edge to an intermediate region inside 
the disk wind, which is < 30 R*, in the studied cases (see Fig . [5j ■ 

According to BP, not all values of k and A are capable of 
producing super- Alfvenic flows at infinity, a condition that en- 
ables the collimation of the disk wind and the formation of a 
jet. The condition necessary for a super-Alfvenic flow is that 
kA(2A - 3) 1/2 > 1, which is true in all plots in this paper. Another 
condition necessary for a solution is that it must be a real one 
and the velocities inside the flux must increase monotonically, 
which more constrains the set of values of a:, A, and ■&, which 
can produce a physical solution. The BP solution is a "cold" 
wind solution, where the thermal effects are neglected and only 
the magneto-centrifugal acceleration is necessary to produce the 
outflow from the disk. The results presented here show that the 
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Fig. 7. Ha profiles when M acc = KT 8 M G yr-' with r magjM Ax = 9000 K and T w 



= 10 000 K. Left panels show (a) the Ha 



profile for different values of A, and (b) the corresponding poloidal velocity profiles for the innermost line of the outflow. The lines 
in the left panels are: A - 30 (solid line), A - 40 (dashed line), A - 50 (dash-dotted line) and A - 60 (dotted line). The right panels 
show (c) the Ha profile for different values of k, and (d) the corresponding poloidal velocity profiles for the innermost streamline of 
the outflow. The lines in the right panel are: k - 0.03 (solid line), k - 0.06 (dashed line), k - 0.09 (dash-dotted line) and k — 0.12 
(dotted line). The mass loss rate is kept constant at 0.1M acc . 



disk wind contribution to the Ha flux dependents very much on 
the disk wind parameters used to calculate its self-similar so- 
lutions. By changing the values of k, A, and -&q, we change the 
self-similar solutions of the problem. Thus, changing the veloc- 
ities and densities inside the disk wind, leads to a displacement 
of the blue-shifted absorption component and a variation, more 
visible around the rest velocity, in the overall disk wind contri- 
bution to the flux. For a "warm" disk wind solution, with strong 
heating mechanisms acting upon the accretion disk surface, and 
producing an enhanced mass loading o n the disk wind, a much 
lower A » 2 — 20 would be necessary (Cas se & Ferreiral l2000). 
We worked only with the "cold" disk wind solution, but it is pos- 
sible to use a "warm" solution, if we find the correct set of values 
for A and k. 



Figure [9] shows a dependence of the line profile with incli- 
nation i similar to that found by KHS using their disk-wind- 
magnetosphere hybrid model. In both models, the Ha line pro- 
file intensity decreases as i increases, and so does the line 
equivalent width (EW). These results contradict the result from 
Appenzel ler et alj ([2005), who have shown that the observed Ha 
from CTTS increase their EW as the inclination angle increases. 
A solution to that problem, as shown by KHS, is to add a stellar 
wind contribution to the models, because that would increase the 
line EW at lo wer inclination. A lar ger sample than the 12 CTTS 
set used bv lAppenzeller et all (120051) should also be investigated 
before more definite conclusions are drawn. However, there is a 
great difference between the blue-shifted absorption position in 
our model and in the KHS model. In their model, the blue-shifted 
absorption position is around the line center, while in ours it is 
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Fig. 8. Ha profiles (a) and poloidal velocity profiles for the in- 
nermost disk wind streamline (b) for different values of the 
launching angle iV The lines are: )?o = 14.6° and a - 0.23 
(gray solid line), § = 33.4° and a - 0.43 (dot-dot-dot-dashed 
line), § = 46.7° and a = 0.63 (dot-dashed line), § = 55.6° 
and a = 0.83 (dashed line). The full black line in (a) repre- 
sents the profile with only the magnetosphere. In all the plots, 
'Mgyr 1 , r mag ,MAX = 9000 K and 
= 10 000 K. 



at a much higher velocity. This difference aris es from the dif- 
ferent disk wind geometry used in both models. iKurosawa et al.l 
(2006) used straight lines arising from the same source point 
to model their disk wind streamlines, which results in the inner- 
most streamlines having a higher launching angle and thus lower 
projected velocity along the line of sight in that region when i is 
large. Instead we use parabolic streamlines, and these lines have 
a higher projected velocity in our default case along the line of 
sight, displacing our absorption feature to higher velocities. 

In both models it is common to have the magnetospheric 
contribution as the major contribution to the Ha flux, with a 
much smaller disk wind contribution. In a few of the profiles cal- 



culated by KHS, the disk wind contribution is far more intense 
than the magnetospheric contribution. With their disk-wind- 
magnetosphere hybrid model and using M acc = 10~ 7 Moyr _1 , 
r wind = 9000 K and p = 2.0, KHS have obtained a disk wind 
contribution to the Ha flux one order of magnitude higher than 
the magnetospheric contribution. In their model, p is the wind 
acceleration parameter, and the higher its value, the slower is the 
wind accelerated, which leads to a denser and more slo wly accel- 
erated outflow in some regions of the disk wind. Kurosawa et al.l 
(2006) used an isothermal disk wind. A large acceleration pa- 
rameter (J3 - 2.0), leading to a large region with a high density, 
which is coupled with an isothermal disk wind at high temper- 
ature (r w i n d = 9000 K), would then lead to a very intense disk 
wind contribution. However, the p = 2.0 value of their model 
is an extreme condition. For our model to produce a disk wind 
contribution as high as in that extreme case, we would also need 
very extreme conditions, although due to the differences between 
the models, these conditions should present themselves as a dif- 
ferent set of parameters. The modified BP disk wind formulation 
used in our model produces a very different disk wind solution 
than that used by KHS, the disk wind acceleration in our solution 
is dependent on A, k and §q parameters (see Figs.[7}3, [7Jl and|8). 
We also use a very different temperature law (see Fig. [2), which 
is not isothermal, and a different geometric configuration. It is 
however also possible to obtain a very intense disk wind contri- 
bution with the correct set of parameters, which corresponds as 
in KHS to a very dense and hot disk wind solution. 

Comparing our resul t s with the classification scheme pro- 
posed by Reipurth et al. (1996), we can see for the Ha emis- 
sion lines of TTS and Herbig Ae/Be stars that most of the pro- 
files in this paper can be classified as Type I (profiles symmetric 
around their line center) or Type III-B (profiles with a blue sec- 
ondary peak that is less than half the strength of the primary 
peak). These two types of Ha p rofiles correspond to a lmost 60% 
of the observed profiles in the Reipurt h et al.l ([1996) sample of 
43 CTTS. Figure @t a ls° shows a profile that can be classified 
as Type IV-R. We have not tried to reproduce all types of pro- 



files in the classification scheme of Reipurth et 



1 types ot ] 
alj JT996), 



but 



with the correct set o parameters in our model, it should be pos- 
sible to reproduce them. Even with our different assumptions, 
our model and KHS model also produce many similar results, 
and both models can reproduce many of the observed types of 
Ha line profiles of CTTS, but the conditions necessary for both 
models to reproduce each type of profile should be different. 

We are working with a very simple model for the hydrogen 
atom, with only the three first levels plus continuum. A better 
atomic model, which takes into consideration more levels of the 
hydrogen atom, would give us more realistic line profiles and 
better constraints for the formation of the Ha line, and would 
allow us to make similar studies using other hydrogen lines. The 
assumption that the hydrogen ground level is in LTE due to the 
very high optical depth to the Lyman lines is a very good ap- 
proximation inside the accretion funnels, and even in the densest 
parts of the disk wind, but it may not be applicable in the more 
rarified regions of the disk wind. In those regions, the interven- 
ing medium will be more transparent to the Lyman lines, thus 
making the LTE assumption for the ground level invalid. Using 
a more realistic atom in non-LTE would lower the population of 
the ground level, and in turn the populations of levels 2 and 3 
would increase. The line source function [Eq. d35l >l depends on 
the ratio between the populations of the upper and lower levels, 
so, if the population of the level 3 increases more than the popu- 
lation of the level 2, the Ha emission will increase, otherwise the 
line opacity will increase. In any case, this will increase the disk 
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wind effects upon the Ha line, thus making the disk wind con- 
tribution in cases with low to very low accretion rates important 
to the overall profile. The real contribution will depend on each 
case and is hard to predict. Still, even with all the assumptions 
used in our simple model, we could study how each of the disk 
wind parameters affects the formation of the Ha line on CTTS. 
Moreover, the general behavior we found might still hold true 
for a better atomic model, only giving us different transition and 
limiting values. 

It is necessary to better understand the temperature struc- 
tures inside the magnetosphere, and inside the disk wind, if we 
wish to produce more precise line profiles. A study ab out the 
temperature structure inside the disk wind done bv lSafierl £1993) 
found that the temperature inside the disk wind starts very low 
(< 1000 K), for lines starting in the inner 1 AU of the disk, then 
reaches a maximum temperature on the order of 10 4 K, and re- 
mains isothermal on scales of ~ 100 - 1000 AU. That temper- 
ature law is very different from the one used in our models, 
and it would produce very different profiles than the ones we 
have calculated here, because the denser parts of the outflow, 
in ISafierl (1 19931). have a much lower temperature than in our 
model. But Safier (1993) only considered atomic gas when cal- 
culating the ionization fraction, and did not consider H2 destruc- 
tion by stellar FUV photons, coronal X-rays, or endothermic re- 
actions with O and OH. A new thermo-chemical study about 
the disk wind in Class , Clas s I, and Class II objects is pre- 
pared by Panoglou et al.| (1201 Oh . where they addr ess thi s prob- 
lem, adding all the missing factors not present in ISafierl {1993) 
and also including an extended network of 134 chemical species. 
Their results produce a temperature profile that is similar to the 
ones we used here, but which is scaled to lower temperatures. 
Nevertheless, the innermost streamline they show in their plots 
is at 0.29 AU, which is around 30 R* in our model and their tem- 
perature seems to increase as it nears the star. This strengthens 
the conviction that the temperature structure used in our model 
is still a good approximation. 

The parameter that defines the temperature structure inside 
the magnetosphere and inside the disk wind is the A parameter, 
defined by Har tmann et al . ( 1982), and it represents the radiative 
loss function. It is then possible to use A to infer the radiative 
energy flux that leaves the disk wind in our models, and even 
in the most extreme scenarios, the radiative losses in the disk 
wind are much lower than the energy generated by the accretion 
process. Our model is consequently not violating any energetic 
constraints. 

In a following paper we plan to use our model to reproduce 
the observed Ho- lines of a set of classical T Tauri stars exhibit- 
ing different characteristics. After that, we plan to improve the 
atomic model to obtain more realistic line profiles of a large va- 
riety of hydrogen lines. 



6. Conclusions 

We developed a model with an axisymmetric dipolar mag- 
netic field in the magnetosphere that was firstly proposed by 
Hartmann et al. ( 1994), and then improved by Muzerolle et al. 
( 1998, 2001), adding a self-similar dis k wind component as pro- 
posed by Blandford & Pavnd dl982l) . and a simple hydrogen 
atom with only three levels plus continuum. We have investi- 
gated which are the most important factors that influence the Ha 
line in a CTTS, and which region in the CTTS environment is the 
most important for the formation of this line. We have reached 
the following conclusions: 
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Fig. 9. Ha profiles with different star-disk system inclinations, 
i, M acc = lO^Moyr- 1 , r mag , M Ax = 8000 K and r wind , MA x = 
10 000 K. For each i, there is a plot showing the total profile 
with both magnetosphere and disk wind contributions (black 
line), and a profile with only the magnetosphere contribution 
(gray line). The lines are: i - 15° (solid), i - 30° (dotted), 
i = 45° (dashed), i = 60° (dot-dashed), and i = 75° (dot-dot- 
dot-dashed). 



- The bulk of the Ho- line flux comes from the magnetosphere 
for CTTS with moderate (M acc « 10~ 8 M Q yr~ l ) to low mass 
accretion rates (M acc < 10 _9 M o yr _1 ). For CTTS with high 
mass accretion rates (M acc > lO^ 1 M Q yr~ l ), the disk wind 
contribution is around the same order of magnitude as the 
magnetospheric contribution. In very extreme scenarios, it 
is possible to make the disk wind contribution much more 
important than the magnetosphere's contribution. 

- There is a transition value of ^j^max below which no blue- 
shifted absorption feature can be seen in the Ha line pro- 
file. This transition temperature depends on the mass loss 
rate, and is lower for higher values of M] oss . This transi- 
tion value can reach from r w j n d,MAX ^ 6000 K for CTTS 
with very high mass loss rates (Mi oss > lO _8 M yr _1 ), to 
?wind,MAX ~ 9000 K for stars with moderate mass loss rates 
(Mioss ~ 10 _9 M o yr _1 ). This transition temperature must 
not be much higher than 10000K, otherwise most of the 
hydrogen will be ionized. The disk wind in CTTS with 
Mi oss < lO^Moyr -1 should not produce any blue-shifted 
absorption in Ha. The other disk wind parameters that affect 
the outflow density can also change the value of this transi- 
tion r w ; n( i,MAX- 

- Most of the disk wind contribution to the Ha flux comes 
from the inner regions of the disk wind, and the blue- 
shifted absorption feature that can be seen in most of the 
observed profiles is also formed around the same region. For 
Mioss ^ 10~ 8 M yr _1 , most of the disk wind Ha flux comes 
from the outflow launched from r < 20 R* in the disk. This 
region is larger for higher mass loss rates. The size of this re- 
gion depends on the mass loss rate, density, and temperatures 
inside the disk wind. 

- The blue-shifted absorption component dependents very 
much on the size of the disk wind, but except for this feature, 
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the disk wind contribution to the Ha line barely changes if 
the Mioss is kept constant while r^ changes. In a scenario 
with a small r& , the disk wind density might surpass a cer- 
tain limiting value (po > 10~ 10 gcrrT 3 ), making the disk 
wind contribution to the profile very strong, even the main 
contribution to the profile in more extreme cases. This value 
for po should change if we instead use a more realistic hy- 
drogen atom. 
- The disk wind parameters that define its angular momen- 
tum (A), mass to magnetic flux ratio (k), and launching an- 
gle (#o) are the main factors that define the position of the 
blue-shifted absorption feature in the line profile. By chang- 
ing these parameters, it is possible to displace this absorption 
feature to higher or lower velocities in the blue-ward part of 
the profile. 

Summarizing, our results indicate that it is necessary to have 
a substantial outflow (Mi oss > 10~ 9 M Q yr~ l ) from the inner 
disk (r < 0.2 AU) for the appearance of a visible overall disk 
wind contribution to the Ha line profile, including the blue- 
shifted absorption feature. For cases with extreme mass loss 
rates (Mi oss > 10~ 8 M & yr~ l ), this region in the inner disk, which 
is the most important for the Ha line formation, should reach 
r > 0.2 AU. The maximum temperature inside the disk wind 
should also be higher than a limiting value, which depends on 
Mi oss , for the blue-shifted absorption feature to appear in the Ha 
line profile. 

We could find better constraints to the model, which would 
solve the degeneracy encountered in our results, by simultane- 
ously fitting other observed hydrogen emission lines of these ob- 
jects. To do this it is necessary to use a better atomic model. This 
model can also be useful to better understand the Herbig Ae/Be 
stars environment, which, like the CTTS, also produce a diverse 
amount of emission l ine profiles (Finkenzeller & M undllll984t 
iMuzerolle et ai1l2004l) . 
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